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Abstract 

Large amplitude kinetic Alfvenon (exact Alfven soliton) matching condition is investigated in 
quasineutral electron-ion and electron-positron-ion plasmas immersed in a uniform magnetic field. 
Using the standard pseudopotential method, the magnetohydrodynamics (MHD) equations are 
exactly solved and a global allowed matching condition for propagation of kinetic solitary waves is 
derived. It is remarked that, depending on the plasma parameters, the kinetic solitons can be sub- 
or super- Alfvenic, in general. It is further revealed that, either upper or lower soliton speed-limit 
is independent of fractional plasma parameters. Furthermore, the soliton propagation angle with 
respect to that of the uniform magnetic field is found to play a fundamental role in controlling the 
soliton matching speed-range. 

PACS numbers: 52.30.Ex, 52.35.-g, 52.35.Fp, 52.35.Mw 
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I. BACKGROUND 



In the collective processes in plasma if the magnetic field perturbations are also taken into 
account, in addition to the ion-acoustic mode, the Alfven mode will also be present. Hence, 
the plasma will be dispersive if the kinetic effects like ion-drift velocity are also included 
in the ion continuity and the ion momentum equations. In fact, it is because of inclusion 
of these effects that the propagating waves are named the kinetic waves. In other words, 
when the perpendicular wavelength is comparable to the ion gyro-radius, the ions will no 
longer follow the magnetic field lines of force, while the electrons/positrons, due to their 
small Larmor-radius, will be attached to the magnetic force lines. This causes an effective 
charge separation which as a result will lead to the kinetic Alfven waves. Therefore, the ion 
drift-velocity is mainly responsible for the Alfven waves kinetic effects. On the other hand, 
kinetic Alfven solitary waves occur due to the interplay between the dispersive behavior and 
nonlinear steepening effects 

Nonlinear Alfven waves such as Alfvenons (exact Alfven soliton)and shocks play impor- 



tant role in solar and astrophysical plasmas |2|, |3(. Stasiewicz [4| has recently suggested that 
Alfvenons can provide an explanation for various aspects of electromagnetic energy dissipa- 
tion and heating in the solar corona and in planetary magnetospheres. Early investigations 
of Alfvenons date back to 1976 . Application of exact methods such as the Sagdeev 
approach is a rigorous tool in the study of the nonlinear phenomena 7H9J in plasmas, which 
provides vital information concerning the wave dynamics. There has been many recent stud- 
ies, mostly using the Sagdeev method, on the possibility of solitary Alfven waves in diverse 
plasmas with dust and positron contaminant, considering the plasma inhomogeneities and 
various charge distribution effects both in low-beta and inertial limits 1CH26I]. 



Investigations show that, both rarefactive and compressive solitary waves exist depending 
on the fractional plasma parameters and the value of plasma-beta (/3 = Bq)i 
a parameter which in most cases is assumed to be smaller than unity. Furthermore the 
theoretical studies of kinetic Alfvenons reveals the existence of both sub- and super- Alfvenic 
nonlinear structures. In a recent study O. P. Sah 27J] has reported the existence of super- 
Alfvenic double-layers in electron-positron-ion plasmas. Woo etal. 28( have found that 
the collisional effects can be important in converting the soliton-shape structures to double- 
layers in dusty plasmas. Prasanta Chatterjee etal. [291 ]. have recently studied the effect 
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of ion-temperature on large-amplitude kinetic solitary waves and concluded that the ion 
temperature has important effect on shape of solitary kinetic Alfven waves and double- 
layers. Current investigation is an attempt to study the most general cases of electron-ion 
and electron-positron- ion three dimensional plasmas consisting of warm electrons/positron 
and ions considering the two dimensional drifts of ions. We will report features in arbitrary 
amplitude soliton dynamics which is common in all electron-ion and electron-positron-ion 
plasmas regardless of the charge distribution function for plasma species. The presentation 
of the article is as follows. In Sec. [Til the basic fluid formalism is presented. The general 
pseudo-potential solution and the required conditions on the solitary propagations is given 
in Sec. Ill II We extend our methodology to the electron-positron cases in Sec. [TV] and give 
several examples in Sec. [V] 



II. MAGNETOHYDRO DYNAMICS MODEL 



We use the conventional magnetohydrodynamics (MHD) formalism for a three- 
dimensional homogenous and collisionless quasineutral plasma present in a uniform mag- 
netic field, B = B e z , which is directed along the z-axis. In this formalism two potentials, 



namely, <p and ip, describe parallel and perpendicular perturbations in the electric field [30j, 
E, respectively. Furthermore, we consider only kinetic perturbations in the magnetic field. 
We also assume that the ions possess both polarization drift and E x B velocity components 
in x-y plane. Hence, the set of normalized MHD fluid equations read as pjj 

d t n e + <9|| (n e u\\ e ) = 0, 

dm + V ± • (njU_u) + <9|| (niU\\i) = 0, 

u±i = E ± x S|| + dtE±, E± = -Vr0, V± = (d x , d y ) , 

Q [d t u\\ e + U|| e <9||U|| e ] - d\\<p + n~ l d\\P e {n e ) = 0, 

d t u\\i + (uj_i ■ V_l) u\\i + u\\id\\u\\i + d\\<p + n^d^P^rii) = 0, 

d\\ [A ± (ip - <p)] = (P/2) d t J h Aj. = (dl, &*) , J\\ = mu\\i - n e u\\ e , 

where, Pj and J are density, velocity, mass, pressure of j-th species (j = e, i) 

and the total current density, respectively. Furthermore, the notations || and _L refer to the 
parallel and perpendicular (with respect to direction of the magnetic field) vector compo- 
nents. The parameter (3 = 2{i n i0 k B T e /BQ (n i0 being the equilibrium ion number density) 
is the ratio of thermal to magnetic pressure and the quantity Q is the electron to ion mass 



ratio to be neglected in our analysis. In obtaining the normalized equation set (Eqs. (JTJ) 
we have used the general scalings 

{x,y,z} ^ ^-{x,y,z}, t -> ±, nj ^ n i0 %, uj ->■ c s u j; 

where, w C j = eBo/m i} c s = ■sj e/rni are ion cyclotron-frequency and acoustic speed, respec- 
tively. The value of parameter e will be defined based on the electronic charge distribution. 



III. KINETIC SOLITONS IN ELECTRON-ION PLASMAS 

Being interested in stationary wave solutions moving at constant velocity, we reduce 
Eqs. ([1]) to co-moving stationary soliton frame by changing into the new coordinate £ = 
l x x + l y y + l z z — Mt {l"i + l1 + l 2 z = 1), with M = V/Va being the normalized matching soliton- 
speed, where Va = B /y/fi niomi is the Alfven wave-speed. By changing the coordinate, 
making use of the quasineutrality condition, n e = rii, and integrating with appropriate 
boundary conditions ( lim ip = 0, lim uj = 0, lim rij = 1), the reduced set of equations 
become 

IxUxi + lyUyi + l z u zi = M(l - n' 1 ) 
u ze = Zj 1 M(l - n- 1 ) 
u xi = Ml x d^ip - lyd^ip 
Uyi = l x d^ + Mlyd^ip 

u zi = l z M- l n<5>*(n), $*(n) = /" n- l d n P{n)dn, 
l z (l 2 x + l 2 y ) (%V - = nf3M {u ze - u zi ) /2, 

where, P{n) = P e {n) + Pj(n) and $*(n) is the effective plasma potential due to the total 
pressure. The following general differential equation is obtained by combining Eqs. ([3]) 



d 2 ip(n) 1 



f l-n) + - — $*(n) - (ri) + (1 - n~ l 



(4) 



2l 2 z 2 M 2 

Algebraic manipulation and integration with the aforementioned boundary conditions, re- 
sults in the well-known energy integral of the form 
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with the desired pseudo-potential 



U(n) = 
where 

<p(n) = 



[d n <p(n)]~ 
I'll 



d(f(n) 
dn 



/3M 2 
~2Zf 



(n - 1) - ^-$» + gd>*(n) + (n- 1 - 1) 



n l d n P e (n)dn. 



dn. (6) 



(7) 



The possibility of solitary excitation relies on some conditions to satisfy, simultaneously, 
namely 



U(n)\ 



dU{n) 



71=1 



dn 



0. 



d 2 U{n) 



n=l 



dn 2 



< 0. 



(8) 



n=l 



It is further required that for at least one either maximum or minimum nonzero n- value, 
we have U(n m ) = 0, so that for every value of n in the range n m > n > 1 (compressive 
soliton) or n m < n < 1 (rarefactive soliton), U(n) is negative (it is understood that there is 
no root in the range [1, n m ]). In such a condition we can obtain a potential minimum which 
describes the possibility of a solitary wave propagation. The stationary soliton solutions 
corresponding to this pseudo-potential which satisfies the mentioned boundary- conditions, 
read as 

dn 



(9) 



y/-2U(n) 

The conditions for the existence of a solitary propagation stated above require that, first, 
it takes infinitely long pseudo-time (£) for the system to get away from the unstable point 
(n = 1). This statement requires that d n U{n) | n =i= or equivalently d^n |{=_oo= in 
parametric space, as it is also inferred by the shape of a solitary wave. Thereafter, moving 
forward in pseudo-time (£ ) axis, the localized density perturbation reaches a maximum or a 
minimum at n — n m (if it exists) at which the pseudo-speed (d^n) of the analogous particle 
bound in pseudo-potential (U(n)) region of 1 > n > n m (or 1 < n < n m ) reaches zero 
again and it returns back. Note that, in the parametric space, from equation Eq. (Q, it is 
observed that in physical situation U (n) should be negative for solitary (non-periodic) wave 
solution, which is clearly satisfied if d nn U{n) | n =i< and U{n m ^ 1) = 0. Note also that, 
both the requirements U(n) | n =i= and d n U(n) |„=i= follow from the equilibrium state 
assumption at infinite pseudo-time (£ = ±oo) before and after perturbation takes place, 
i.e. d^n |^ = ±oo= d^n \(=±oo= 0. However, there is a special case with d n U{n) \ n=Um = for 



which the density perturbation is stabilized at the maximum or minimum density n = n m 
(the analogous particle never returns back). This situation regards to the existence of 
a double-layer in plasma which is not considered here. Extra requirement such as 
(n m — l)d^n |^=o< given in Ref. jsi| is needed for exclusion of the possibility of 
double layer or shock-like structure in addition to the ones presented in Eq. (151). 
However, in the present analysis we only attempt to give a general criterion 
for the soliton matching speed with the assumption that such density profile 
exists for the set of given plasma parameters in Eq. ((6j). Therefore, we continue 
to evaluate the matching condition with the assumptions that a root n m ^ 1 
exists and (n m — l)d^n |g=o< (to exclude shock-like solutions). Moreover, the 
pseudopotential given by Eq. (jH]) and its first derivative vanish at n — 1, as required by 
the two first conditions in Eq. (jSJ). Also, direct evaluation of the second derivative of the 
Sagdeev potential, Eq. fl6]), at unstable point, n = 1, leads to 



Thus, in order for the existence of oblique kinetic Alfven solitary propagations, the matching 
soliton-speed must satisfy the following general inequalities 



This condition is analogous to the one presented for oblique electrostatic solitary propaga- 
tions in magnetoplasmas given in Ref. [32]. However, in current case the propagation angle 
with respect to the magnetic field, l z , plays a fundamental role in controlling the matching 
speed range. Note that there is also a critical /3-value (f3 cr = 2 [d n $*(l)] _1 ) at which no ki- 
netic Alfvenons can propagate. Furthermore, no kinetic Alfvenons are allowed to propagate 
at right angle to the magnetic field. 

IV. EXTENSION TO ELECTRON-POSITRON-ION PLASMAS 

We extend the methodology used above to obtain matching conditions for electron- 
positron-ion kinetic Alfven solitary propagations. To this end, we change to the new coordi- 
nate £ = l x x + l y y + l z z — Mt {l"i+ly+l 2 z = 1), use appropriate boundary conditions lim = 




(10) 




(11) 
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and lim nj jejP = {l,a,a — 1}, where a = n e o/riio > 1 (the special case of a = 1 leads to 
the quasineutral electron- ion magnetoplasma) , and employ the quasineutrality condition, 
n i = n e — n p , to obtain similar relations as Eqs. (EJ) 



^x«xi + + = M(l - nj' 1 ) 
n e u ze = l~ l M{n e - a) 
n p u zp = l~ l M(n p - a + 1) 
u xi = Ml x d^ip - lyd^ip 
Uyi = l x d^ip + Ml y d^i> 

u zi = l z M- l ni^*{ni), ^*(rii) = <p(in) + J™ 1 nj 1 d^P^n^d^ 
h (ll + l 2 y ) (dtfip - dtfcp) = (3M (n e u ze - n p u zp - niU zi ) /2, 



(12) 



where, the function tp(rii) is given by the following relations 



dip 
dn e . 



1 dP e (n e 
n P dn P 



dip 
dn n 



1 dP p (n p ) 
n p dn p 



T^e Tl p ~ T^ii 



(13) 



Using the above relations, we expect to find the function <p = ip(jii) which obviously is 
nontrivial in most cases. However, as it will be revealed, we need to calculate the quantity 
d 9 ni(ip), instead, which is feasible in many cases. Therefore, for the current case, Eq. (j4]) 
changes to 



d 2 (p{nj 
dE? 



I — ll 



(3M 2 



II, 



(ru 



l]ni 
M 2 



y*( n .) + ( n . _ i) 



(14) 



The generalized Sagdeev pseudopotential U(rii), then, reads as 



I -11 



Hi 



dipjrij 
dm 



^{ni - 1) - + j^^in,) + (1 - m) 



drii. 



(15) 



The matching conditions to be evaluated is of the similar shape as before 



< M < ^4^(1); /3d ni **(l) > 2 
however, with the new definition of the effective potential given below 



(16) 



n.i=l 



(17) 
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V. CASE STUDY FOR SOLITON MATCHING CONDITION 



It clearly remarked from Eqs, [TT] and [16] that both sub- and super- Alfvenon may exist 
depending on the equation of state for electrons/positrons and ions. Considering the pressure 
of ions to be Piirii) = n]k B Ti/e, where 7 = (/ + 2)// is the adiabatic constant and / is the 
ion degrees of freedom (DoF) and with the cases of 7 = 1 for isothermal-ions and 7 ^ 1 
for adiabatic-ions, we employ different electronic distributions to evaluate the possibility 
condition for solitary propagations. For instance for the Maxwell-Boltzmann, Lorentzian 



(Kappa), Thomas- Fermi and Fermi-Dirac 33] electron/positron distributions (appropriately 
normalized), we have 



Maxwell { \ e = k B T e 5 = ^ 

n p = {a- l)e- s * 



n e = a(l-v/(k-3/2))- k+1/2 
(a-l)(l + 6cp/(k-3/2))- k+l/2 



Lorentz { e = k B T e 5 = ^ 



n e = a(l + <f) 



.3/2 



n p 



(a - 1)(1 - 5<ff 2 



Fermi ( ' _ ._ \ e = k B T Fe 5 = ^ E 



if = y/l + a 2 /% e 2 / 3 ^ -y/l + a^rg 



Dirac I ' * / — * / \ e = m e c 2 t]l = jf- 



No 



(18) 



where, a is the electron-to-ion equilibrium number- density ratio and Tp e ,p is the elec- 
tron/positron Fermi-temperature, respectively. The parameter t]q is related to the plasma 
mass-density (of white dwarf, for instance, where the relativistic degeneracy occurs) through 
the relation p ~ 2m p no or p(gr/cm 3 ) = po^o with po(gr / cm 3 ) ~ 1.97 x 10 6 , where, m p is the 
proton mass and the cases p{= p/po) "C 1 and p 1 correspond to the nonrelativistic and 



ultrarelativistic degeneracy limits, respectively [33] . The density p is exactly in the range 
of a mass-density of a typical white dwarf (the density of typical white dwarfs can be in 
the range 10 5 < p(gr/cm 3 ) < 10 9 ). The normalized matching Alfvenon-speed in the case of 
electron-positron- ion magnetoplasma is bounded through inequality of the form, Eq. (I16p . 
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with the values of d ni ^f*(l) given as 
Maxwell d ni tf*(l) = 7<r + [a + 5 (a - o 



21 



Lorentz = 7^ + §Ef [« + <K<* - 
Fermi <*rn**(l) = 7^ + |[« + 5(a - a = TF e 

Dirac = 1 ■ a ~ 

V3(a-l) 1/3 V /l +(«- 1 ) 4/3 '?o+ 3Ql/3 V /l + a4/3, ?o 



(19) 



where, the a = 1 case corresponds to the electron-ion limit. For instance, the a = 1 
(electron-ion) limiting case of the Kappa distribution in Eq. (|T9"|) has been investigated 
in Ref. J26]. Note also that for the limit of spectral index k — > 00 the result for the 
Maxwell-Boltzmann distribution is obtained. Evaluation of Eq. (|19j) indicates that the 
effects of adiabatic index, 7, and fractional ion-temperature, a, is to shrink (widen) the 
solitary matching speed-range, while, the effects of fractional electron number-density, a, 
and fractional positron temperature, 5, is to widen (shrink) the corresponding range for 
0d n .^*{l) < 2 (/3d n ^*(l) > 2). It is also observed from the values given in Eq. flTHJ that, 
the increase in the plasma mass-density in a Fermi-Dirac relativistically degenerate plasma 
shrinks (widens) the Mach-number range for the case (3d ni ^*(l) < 2 (j3d ni ^*(l) > 2). 



9 



[1] Roychoudhury, R. and Chatterjee, P., Phys. Plasmas, 5 3828(1998). 

[2] P. Louarn, J. E. Wahlund, T. Chust, HdeFeraudy, A. Roux, B. Holback, P. O. Dovner, A. I. 
Eriksson and G. Holmgren, Geophys. Res. Lett. 21 1847(1994). 

[3] P. O. Dovner and G. Holmgreen, Geophys. Res. Lett. 21, 1827(1994). 

[4] K. Stasiewicz, Phys. Rev. Lett., 96 175003(2006). 

[5] Akira Hasegawa and Kunioki Mima, Phys. Rev. Lett., 37 690(1976). 

[6] M. Y. Yu and P. K. Shukla, Phys. Fluids, 21 1457(1978). 

[7] A. A. Vedenov, E. P. Velikhov, and R. Z. Sagdeev, Yad. Sintez 1 82(1961). 

[8] R. C. Davidson, Methods in Nonlinear Plasma Theory, Academic Press, New York, (1972). 

[9] R. Z. Sagdeev, 1966 "Reviews of Plasma Physics", Vol. 4 ed. M. A. Leontovich (New York: 
Consultants Bureau). 

[10] K.P.Das and F. Verheest,J. Plasma Phys. 41, 139(1989). 

[11] De-Jin Wu and De-Yu Wang, Phys. Plasmas, 2 4476(1995). 

[12] H. Kakati and K. S. Goswami, Phys. Plasmas, 5 4229(1998). 

[13] Xu-yu Wang, Xue-yi Wang and Zhen-xing Liu, Zhong-yuan Li, Phys. Plasmas, 5 4395(1998). 

[14] Anup Bandyopadhyay, Phys. Plasmas, 7 3227(2000). 

[15] Yao Chen, Zhong-yuan Li, Wei Liu and Zhi-Dong Shi, Phys. Plasmas, 7 371(2000). 

[16] RAJKUMARROY CHOUDHURY, J. Plasma Phys., 67 199(2002). 

[17] M. Ansar Mahmood, Arshad M. Mirza, P. H. Sakanaka and G. Murtaza, Phys. Plasmas, 9 
3794(2002). 

[18] H. Saleem and S. Mahmood, Phys. Plasmas, 10 2612(2003). 

[19] M. Arshad Mirza, M. Ansar Mahmood and G. Murtaza, New J. Phys., 5 116(2003). 

[20] L. Yang and D. J. Wu, Phys. Plasmas, 12 112901(2005). 

[21] Cheong Rim Choi and D.-Y. Lee, Phys. Plasmas, 14 052304(2007). 

[22] Yang-fang Li, D. J. Wu and G. E. Morrill, Phys. Plasmas, 15 083703(2008). 

[23] Tao Xu, Bo Tian, Li-Li Li, Xing L and Cheng Zhang, Phys. Plasmas, 15 102307(2008). 

[24] N. Rubab, N. V. Erkaev and H. K. Biernat, Phys. Plasmas, 16 103704(2009). 

[25] W. Masood, H. A. Shah, N. L. Tsintsadze and M. N. S. Qureshi, Eur. Phys. J. D, 59 413(2010). 

[26] Runmoni Gogoi and Manoranjan Khan, Phys. Plasmas, 17 112311(2010). 



10 



[27] O. P. Sah, Phys. Plasmas, 17 032306(2010). 

[28] M. H. Woo, C.-M. Ryu and C.-R. Choi, Phys. Plasmas, 17 053707(2010). 

[29] Prasanta Chatterjee,Taraknath Saha, S. V. Muniandy and C. S. Wong, Phys. Plasmas, 16 
103702(2009). 

[30] B. B. Kadomtsev "Plasma Turbulence", Academic New York (1965) p.82. 

[31] De-Jin Wu, Guang-Li Huang and De-Yu Wang, Phys. Plasmas, 8 2879(1996). 

[32] M. Akbari-Moghanjoughi, Phys. Plasmas, 18 032103(2011). 

[33] M. Akbari-Moghanjoughi, Phys. Plasmas, 18 012701(2011). 



11 



